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ABSTRACT 

Lyutikov (astro-ph/0503505) raised a valid point that for shock deceleration of a highly mag- 
netized outflow, the fate of the magnetic fields after shock crossing should be considered. However, 
his comment that the deceleration radius should be defined by the total energy rather than by 
the baryonic kinetic energy is incorrect. As strictly derived from the shock jump conditions in 
Zhang & Kobayashi (2005), during the reverse shock crossing process the magnetic energy is not 
tapped. As a result, the fireball deceleration radius is defined by the baryonic energy only. The 
magnetic energy is expected to be transferred to the circumburst medium after the reverse shock 
disappears. The evolution of the system then mimicks a continuously-fed fireball. As a result, 
Lyutikov's naive conclusion that the forward shock dynamics is independent on the ejecta content 
is also incorrect. The shock deceleration dynamics and the reverse shock calculation presented 
in Zhang & Kobayashi (2005) are robust and correct. 

In our recent paper Zhang & Kobayashi (2005, hereafter ZK05), we calculated the gamma-ray burst 
(GRB) reverse shock emission when a highly relativistic outflow with a wide range of magnetization parameter 
a interacts with a constant density medium (ISM). We found that throughout the period when the reverse 
shock crosses the ejecta, the lab-frame Poynting flux energy is essentially unchanged. Considering a fireball 
with a total energy E — Ek,q + -Ep.o = Ek.q(1 + where Ek.o — -foMoc 2 is the initial baryonic kinetic 
energy in the ejecta and Ep$ is the initial Poynting flux energy, we drew the conclusion that the radius at 
which the fireball collects 1/To of the initial fireball rest mass, i.e. 

R-y ~ Rdec ~ ( 2 ~ 9 ) ' (1) 
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defines the radius for fireball deceleration in the "thin shell" regime (Eq.[41] in ZK05) when a reverse shock 
exists. This is smaller than the conventional deceleration radius (defined by E, in the a — limit) by a 
factor (1 + 0-) 1 / 3 . As a result, the forward shock bolometric emission level right after this deceleration radius 
is correspondingly fainter by a factor of (I + a) comparing to a — case. 

In a recent astro-ph comment by Lyutikov (2005) , he raised a valid point that the fate of magnetic fields 
after shock crossing, which has been ignored in the version 2 of our paper posted in astro-ph, should be 
considered. The Poynting energy (i?p,o) should be eventually transferred to the ISM. We have incorporated 
this comment in our final version to appear in ApJ (version 3 in astro-ph), and discussed the possibility of 
transferring -Ep.o to the ISM after the shock crossing stage. 

However, there exist severe errors in Lyutikov's comments, which would mislead the readers if not 
corrected. 

The major error in Lyutikov (2005) is that he defines deceleration using the total energy E regardless 
of the a value. This gives 

C ~ (^r^) = ^ec(l + 0-y/ 3 . (2) 
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Such a definition is purely based on his intuition without investigating the rigorous shock jump conditions 
carefully. The consequence of adopting such an assumption unfortunately leads to self-contradictory results 
in the high-er regime. We now elaborate it in the following. 

There are several characteristic radii when discussing the deceleration of a fireball by ISM through 
reverse shocks (Sari & Piran 1995). Besides R 1 , the relevant ones also include Ra, the radius at which 
the reverse shock crosses the ejecta shell, and Rn, the radius at which the reverse shock upstream and 
downstream becomes relativistic to each other. When a increases, a reverse shock forms when the forward 
shock pressure exceeds the magnetic pressure in the outflow (Eqs.[31] and [43] in ZK05). The critical radii 
discussed above may deviate from the values in the a = limit. A detailed analysis shows that for a fixed 
total energy E, Rn does not depend on a, and Ra is roughly smaller by a factor c 1 / 2 (Eqs [36] and [38] 
in ZK05). The relations among these critical radii are characterized by Eqs. (46) or (51) in ZK05. One 
can then categorize various parameter regimes in the T/t 7 — a plane as long as the reverse shock forming 
condition is satisfied (see Fig. 4 of ZK05), where T is the duration of the burst, t 7 is a typical time scale 
defined by the total energy E and the initial Lorentz factor 70 (Eq.[48] in ZK05, where 70 was denoted as 
74)- 

Regardless of the a value, as long as T is long enough (Region (I) in ZK05), one should be in the regime 
of Rn < i? 7 < Ra, i.e. the so-called thick shell regime (cf. Sari & Piran 1995 for the case in the a = limit). 
In this regime, the reverse shock upstream and downstream becomes relativistic with each other before the 
reverse shock crosses the shell. As a result, the ejecta is significanly decelerated after one shock crossing. In 
this regime Ra defines the radius of fireball deceleration. R 1 is no longer meaningful. 

As T becomes smaller, one enters the thin shell regime. In this regime, shock crossing only mildly 
decelerates the shell, and the shell starts to decelerate significantly when it feels a large enough inertia from 
the ISM (at Rdec)- In this regime Rd ec < Rn should be satisfied. In the a — limit, this naturally happens 
when T < t 7 . When a higher a value is considered, one still gets a consistent picture when Rdec is defined 
as R 1 (Eq.[l]), and the transition naturally happens at T < t 1 Q(a), where Q(a) is a function of a (oc a 2 / 3, 
at large a). However, if one uses the un-justified deceleration radius defined by Lyutikov, R^ ec (Eq.[2]), one 
reaches an unphysical regime Rn < Ra < Rdec when o -2 ^ 3 <T/t~ ( < a 2 is satisificd. In such a regime, the 
reverse shock upstream and downstream becomes relativistic first (at Rn), so that after the reverse shock 
crosses the shell (at Ra), the shell should be significantly decelerated. However, the fireball has not reached 
Rdec Y et - the radius Lyutikov expects the fireball to decelerate. This conclusion is obviously absurd, which 
already invalids his definition. Below we will explicitly analyze the inconsistency in Lyutikov's argument. 

Defining R 1 as Rdec in the thin shell regime by ZK05 is not an assumption. It stems from solid energy 
conservation physics and rigorous shock jump conditions. Let us analyze the deceleration process in detail. 

(1) Before the shock crossing process is over, the system contains four regions separated by two shocks 
and one contact continuity. One can write the energy conservation during the process. Assuming that the 
ISM mass collected by the fireball at the end of shock crossing is Mism, The total energy of the whole system 
before the shock crossing is joMqc 2 + Ep^ + Mismc 2 . After shock crossing, the total energy of the system 
is 7(734Mo + 7-Mism)c 2 + Ep, where 7 is the bulk Lorentz factor at the end of shock crossing, 734 is the 
reverse shock Lorentz factor which is ~ 1 in the thin shell case we are discussing, the term 7-Mismc 2 takes 
into account the rest ISM mass and its thermal energy (which is (7 — 1)Mismc 2 ), and Ep is the lab-frame 
Poynting energy after the shock crossing. The energy conservation for the whole system then reads 

(70 - l)M c 2 + {E Pfi - Ep) = ( 7 2 - l)M ISM c 2 . (3) 

According to Eq.(40) of ZK05, strict shock jump condition analysis gives Ep fi ~ Ep in the high-cr regime. 
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The physical reason is that during the shock crossing the magnetic fields are compressed so that their 
comoving energy increases. In the meantime, the whole system is decelerated, so that the lab frame Poynting 
energy Ep = jU' B (where U' B is the comoving magnetic energy) remains essentially unchanged. This means 
that the term (Ep t o — Ep) drops out from the Eq.(3). The remaining equation is the standard fireball 
deceleration equation with an effective fireball energy Ek.o- The fireball starts to decelerate when Mism ~ 
Mo/ '70 is satisfied, at the radius defined by i? 7 (eq.[l]), not at i?^ ec defined by Lyutikov (eq.[2]). 

(2) What is the status of the ejecta beyond i? 7 ? What happens after the reverse shock disappears? Is 
R-dec meaningful at all? 

Again one can write down the energy conservation of the whole system during the process. Let's simply 
re-define 70, M , Ep, as the bulk Lorentz factor, effective mass and Poynting energy at R 7 (Eq.[l]), and 
7, Ep are the Lorentz factor and Poynting energy at a later radius R. We also re-define Mism as the ISM 
mass collected during the process (from R 7 to R). Then the energy conservation can be still expressed as 
equation (3). We show below that 7 must be significantly smaller than 70. In other words, the fireball is 
decelerated beyond R~ r 

If multi-crossing of the reverse shock still happens beyond R 1 , the (Ep t o — Ep) term still drops out from 
the energy conservation equation, so that the fireball must be decelerating. The energy transfer process only 
happens when the reverse shock disappears. 

During the energy transfer epoch (after the reverse shock disappears), fireball deceleration is still going 
on. This is manifested by the argument below. First, it is obvious that one should always have (Ep y0 — Ep) > 
0, i.e. the Poynting energy should only decrease and gets transferred to the medium. Next, the lab frame 
Poynting energy could be written as 

,B' 2 R 3 B' 2 ^B' 2 

Ep = ^ = ^ v ^ = ^T^ = R ^' (4) 

where V is the comoving volumn, which can be expressed as R 3 /j in the spreading phase (which is the case 
in the Regime III of Fig. 4 of ZK05). For the non-spreading case, a similar conclusion could be also drawn. 

If the GRB outflow does not decelerate from R 7 to R^/ ec , as Lyutikov hopes, one would maintain a 
forward shock with a constant thermal pressure in the shocked region. The magnetic pressure B' 2 /8tt 
behind the contact discontinuity should keep constant in order to maintain hydrodynamical equilibrium. 
Eq.(4) then leads to an absurd conclusion, i.e. the lab- frame Poynting energy increases with time (by a huge 
factor). This again invalidates the assumption of Lyutikov that the outflow decelerates at R% ec - 

Then how is the Poynting energy transferred to the ISM eventually? This is an interesting problem 
and needs further detailed investigation. The treatment in ZK05 (ignoring this effect) is a reasonable 
approximation shortly after the reverse shock peak, but will become less rigorous when the reverse shock 
disappears because of the energy transfer from the Poynting energy to the ISM. According to the energy 
conservation equation (3), one can sec that the dynamics depends on how the term (Ep_o — Ep) evolves with 
7. In any case, the system would mimick a continuously-fed fireball, with the Poynting energy gradually 
injected into the system. Physically, B' 2 /8n is related to the pressure at the contact discontinuity during the 
deceleration phase. Right after the shock crossing (at or beyond R y for thin shells and at Ra for thick shells), 
the pressure is balanced at the contact discontinuity, and there is no net energy transfer from the Poynting 
energy to the ISM kinetic energy. As the system decelerates, the magnetic pressure behind the contact 
discontinuity becomes stronger than the thermal pressure in front, so it does work to the shocked region in 
front of the contact discontinuity. This makes the system decelerates less significantly than it would have 
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been without such a magnetic pushing. Such a process keeps going on, until most of the Poynting energy 
is transferred to the ISM. The late afterglow of a high-cr flow could be comparable to that of a low-cr flow, 
but in the early afterglow phase, the afterglow level is significantly lower. The crucial point here is that 
energy transfer phase is separated from the reverse shock deceleration phase, so that the bolometric reverse 
shock emission peak epoch is separated from the bolometric forward shock emission peak epoch. Detailed 
deceleration process is under investigation. One thing is clear: nothing special happens at the radius R^ ec . 

It is now evident that another statement in Lyutikov (2005), i.e. "energy in the forward shock is 
determined by the total energy and is virtually independent of the content of the ejecta" , is also incorrect. 
According to the above analysis, the higher the a, the more prominant the injection process would be. 
Different a leads to different deceleration dynamics and different afterglow lightcurves. The "universal" 
dynamics 7 oc i?~ 3 / 2 suggested by Lyutikov (notice he used oc t~ 3 / 2 so that his t is not the observer's time) 
is only valid when the energy transfer process is over, i.e. when essentially all the energy of the system is 
given to the ISM. However, the energy transfer is a complicated process, and it obviously varies when a 
varies. The naive treatment of Lyutikov (2005, and his previous work) lacks solid physical justification. 

In conclusion, the dynamics of shock deceleration of a magnetized, relativistic outflow presented in ZK05 
is robust. The calculations of the reverse shock emission (which is the main subject of our paper as reflected 
from the title) arc correct. The forward shock emission after the reverse shock disappears is not treated in 
detail, and will be studied more carefully in a future work. We thank Lyutikov 1 for pointing out the flaw 
in our previous version of ignoring the fate of the magnetic fields after shock crossing, and we have added 
a relevant discussion in the final version of the paper to appear in ApJ (version 3 in astro-ph). Yet, the 
flaw does not influence the bulk of the rigorous calculations presented in the paper. In contrast, the main 
arguments presented in Lyutikov (2005) are incorrect, as has been explicitly explained in this reply. 
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1 There have been extensive email discussions between Lyutikov and us about the subject. Later we decided to temporarily 
withdraw the paper from ApJ to add in a discussion about the caveat of the forward shock calculation in our paper, and 
we planned to update the astro-ph version after the paper is finalized. We have also acknowledged Lyutikov for his critical 
comments, and let him know that we are preparing the final version and will send it to him for comments when it is completed. 
It is to our surprise that Lyutikov still decided to expose his incorrect opinion at astro-ph. Although astro-ph might be a chat 
board to exchange ideas, we believe that refereed journals are more appropriate places to lay out scientific view points. We will 
not reply to any further comments on this subject in astro-ph. 
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